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Based on a class of exactly solvable models of interacting bose and fermi liquids, we com¬ 
pute the single-particle propagators of these systems exactly for all wavelengths and energies 
and in any number of spatial dimensions. The field operators are expressed in terms of bose 
fields that correspond to displacements of the condensate in the bose case and displacements 
of the fermi sea in the fermi case. Unlike some of the previous attempts, the present at¬ 
tempt reduces the answer for the spectral function in any dimension in both fermi and bose 
systems to quadratures. It is shown that when only the lowest order sea-displacement terms 
are included, the random phase approximation in its many guises is recovered in the fermi 
case, and Bogoliubov’s theory in the bose case. The momentum distribution is evaluated 
using two different approaches, exact diagonalisation and the equation of motion approach. 

The novelty being of course, the exact computation of single-particle properties including 
short wavelength behaviour. 


I. INTRODUCTION 


Recent years have seen remarkable developments in many-body theory in the form of an assortment of 
techniques that may be loosely termed bosonization. The beginnings of these types of techniques may be 
traced back to the work of TomonagaS and later on by LuttingeiH and by Lieb and Mattisi. The work 
of SawadaE and Arponen and Pajannei in recasting the fermi gas problem in a bose language has to be 
mentioned. Arponen and Pajanne recover corrections to the Random Phase Approximation (RPA) of Bohm 
and Pinesi in a systematic manner. In nuclear physics, bosonization is widely used to study collective 
properties, for an introduction see the book by lachello and Arimai. In the 70’s an attempt was made by 
Luther0 at generalising these ideas to higher dimensions. Closely related to this is the work by Sharp et. 
al.0 in current algebra. More progress was made by HaldaneEl which culminated in the explicit computation 
of the single-particle propagator by Castro-Neto and FradkinEl and by Houghton, Marston et.a 10 and also 
by Kopietz et. al0. Rigorous work by Frohlich and MarchettS is also along similar lines. Also the work 
of Frau et. al§ on algebraic bosonization is relevant to the present article as the authors have considered 
effects beyond the linear dispersion in that article. The exactly solvable models of Calogero and Sutherland 
are of relevance here as well, the exact propagators of these models have been computed by various authors0. 
Recently, these types of models have been generalised to more than one dimension by Ghoslli. 
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The attempt made here is to generalise the concepts of HaldaneEi to accomodate short-wavelength fluc¬ 
tuations where the concept of a linearised bare fermion energy dispersion is no longer valid. To motivate 
progress in this direction, we hnd that it is necessary to introduce two different concepts, one is the canonical 
conjugate of the fermi/bose density distribution, the other is the concept of sea/condensate displacements. 

Histrorically speaking, the idea that the velocity operator could serve as the canonical conjugate of the 
density has been around for a long time, and this has been exploited in the study of Hell by Sunakawa et. 
al.0. However, the authors are not aware of a rigorous study of the meaning of this object, in particular, an 
explicit formula for the canonical conjugate of the density operator has to the best of the authors’ knowledge 
never been written down in terms of the held operators. The work by Sharp et. ahlil comes close to what 
we are attempting here. 

The concept of a sea-displacement is a generalisation of the traditional approach used for bosonizing 
one-dimensional systems such as the Tomonaga-Luttinger0,i models. There, one introduces bose helds that 
correspond physically to displacement of the fermi surface (in ID, fermi points). These bose helds have 
simple forms relating them to number-conserving product of fermi helds. The held operator is obtained by 
exponentiating the commutation rule between the surface-displacement operator and the held operator. By 
analogy, we generalise these ideas, so that one is no longer restricted to be close to the fermi surface. The 
way this is done is to postulate the existence of bose helds that correspond to displacements of the fermi-sea 
rather than just the fermi surface. From this it is possible to write down formulas for the number-conserving 
product of two fermi-helds in terms of the bose helds. A similar construction is possible when the parent 
helds are bosons, but here, we hnd that instead of sea-displacements, we have to introduce operators that 
correspond physically to displacements of the condensate. Actually, the bose case is much simpler and a 
mathematically rigorous formulation of this correspondence is possible. This is a boon, since we use this fact 
and make plausible the analogous correspondence in the fermi case. The assertions in the fermi case are not 
proved ’’rigorously”, rather are made exceedingly plausible by analogy. This is the main drawback of this 
article. 

This article is organised as follows. In the next section, we present some formulas that relate the number 
conserving product of two fermi/bose helds to the relevant sea/condensate-displacement operators that are 
postulated to be canonical bosons. The sea/condensate-displacement operators in turn may be related to 
the parent fermi/bose helds, as it happens, this formula is simple in the case when the parent helds are 
bosons but is difficult in the case when the parent helds are fermions. 

Following this, we write down a generic formula for the fermi/bose held-operator in terms of the density 
operator (operator-valued distribution, to be precise) and its canonical conjugate. The new ingredient in 
this section is the canonical conjugate on the density operator. This quantity may in turn be related to 
currents and densities. We hnd that these formulas are ambiguous unless a proper choice is made for a 
certain phase functional. For bosons, we hnd that this choice is the zero-functional but for fermions it has to 
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be determined by making contact with the free theory (done in the section following the one just described). 

Combining the two previous sections, we write down in the next section, formulas for currents and densities 
in terms of the sea/condensate-displacements, the field-operator has a formula in terms of the sea/condensate- 
displacements as well. Contact is made with the propagator of the free theory and the undetermined phase 
functional of the previous sections is determined for the fermi case. In the next section, interaction terms are 
introduced that correspond to two-body repulsive interactions. It is argued and demonstrated that selectively 
retaining parts of the interaction that are quadratic in the sea/condensate-displacements amounts to doing 
Bogoliubov/RPA theory. Corrections to this quadratic hamiltonian are easy to write down but are not used 
to compute corrections to RPA/Bogoliubov theory as this requires significantly more effort. It is found that 
the diagonalisation of the RPA hamiltonian is rather tricky if one wants to recover both the particle-hole 
modes and the collective mode. In the end, closed formulas are written down for the fermi propagator in all 
three spatial dimensions and their various qualitative features are examined. This completes the solution of 
the many-body problem in the RPA/Bogoliubov-limit. 

The appendices are as follows. Appendix A contains a detailed proof of the correspondence between the 
number-conserving product of two bose fields and the corresponding condensate-displacements. Appendix 
B involves writing down similar ideas for fermi systems. However here, the various assertions are only made 
plausible unlike in the bose case where a rigorous solution is possible. Appendix C is devoted to proving the 
assertion that retaining only terms linear in the sea-displacements in the definition of the density recovers 
the RPA. Appendix D involves a derivation of the formula for the momentum distribution of the ID system 
using the equation of motion approach. Appendix E contains some technical statements regarding the proof 
of uniqueness of the formula relating the fermi field with the corresponding currents and densities. 


II. EXPRESSING PRODUCTS OF PARENT FIELDS IN TERMS OF SEA-DISPLACEMENTS 

In this section we introduce canonical bose fields called sea-displacements in the fermi case and condensate- 
displacements in the bose case. First, we write down a formula for the number conserving product of two 
bose fields in terms of the condensate displacement operators. A rigourous proof of this is relegated to 
Appendix A. The correspondence is made plausible by making several observations about these formulae. 
Let us first focus on the bose case. Let 6 q and 6 )^ be canonical bose operators. From these, we may construct 
other bose operators defined as follows(q 7 ^ 0 ), 

c?q/2(q) = {^=)blbq (1) 

V-'Vo 

and, 

do(0) = 0 (2) 
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where Nq = b^bo. This is the condensate displacement annhilation operator. It is so named for the following 
reason. The definition suggests that this operator removes a boson from among those that are not in the 
condensate and returns it to the condensate, thereby displacing the latter. The reason for the redundant 
momentum label in the notation (iq/2(q) becomes clear if one realises that a more general object would 
be a sea-displacement annhilation operator <ik+q/2(q) • Since the condensate corresponds to k = 0 , we 
have just the condensate displacement annhilation operator. In fact, it will be shown subsequently that 
for the fermi case we have to deal with this more general object namely, the sea-displacement annhilation 
operator. It may be shown that (see Appendix A) this object dq/2(q) satisfies canonical bose commutation 
rules. Also a formula is possible for the number-conserving product of two parent bosons in terms of these 
condensate-displacements. The formula is written down below and proved in Appendix A. 

^lk-(-q/2^k-q/2 = A"ol5k,0<5q,0 + ['5k+q/2,0 ( \/^)<^k( —Q) + ^k-q/2,0<^k (l) (\/^^)] 

+ ^(l/2)(k+q/2)('^ + q/ 2 )c^(l/ 2 )(k-q/ 2 )(k - q/2) 

where. 

No = N 


Mq/2(q), A^] = 0 

( 5 ) 

iV = ^6^6q 

(6) 


q 


also the object do(0) = 0, by definition. 

The way the authors initially derived this formula is as follows. One starts off with the observation 
that the object &k^.q/2^k-q/2 is the only one that enters in the hamiltonian of number-conserving systems. 
Furthermore, it satisfies closed commutation rules amongst other members of its kind. One is therefore led 
to look for formulas for these objects in terms of other bosons with a view to make the full hamiltonian more 
easily diagonalisable. In particular, if there were bose operators c?q/2(q) and dq/2(‘l) snch that ^k+q/2^k-q/2 
was exactly linear in these bosons, then the full hamiltonian would indeed be exactly diagonalisable. We 
find that this is not the case and there are corrections to this linear term and it so happens that introduction 
of a quadratic term in the condensate-displacements in fact makes the correspondence exact. The authors 
are not aware of a deeper reason behind this simple formula that terminates after including the quadratic 
term, after all, the formula for the parent annhilation operator 6q in terms of the condensate displacements 
is formidable as we shall soon see. The bose case being so simple and exact can be used as a benchmark 
to write down corresponding formulas in the fermi case, where rigorous proofs are much harder to come by. 


( 3 ) 

( 4 ) 


4 



The authors also have in mind generalisations to relativistic systems, where one might profit by following 
the above prescription. In particular, it would be fascinating to see if the ideas above were useful in getting 
non-perturbative information regarding gauge theories like QED, QCD e.t.c. But this is far into the future. 
For now, let us try to write down a similar correspondence for the nonrelativistic fermi system. 

As mentioned earlier, for fermi systems, it is necessary to postulate the existence of a sea-displacement 
annhilation operator, denoted by ak(q). A formula for this in terms of the fermi fields is extremely difficult 
to deduce. In Appendix B, attempts are made to do exactly this. There it is pointed out that these 
objects satisfy canonical boson commutation rules. The important issues that enable one to draw practical 
conclusions, fortunately do not depend very much on the technical details. In Appendix B and in the sections 
that follow, we show how to extract the necessary physics while cirumventing the technical details. It must 
be pointed out however that this drawback is regrettable. Let us merely quote the final answers and later 
on make these formulas plausible. The RPA form of the number conserving product of two fermi fields in 
terms of the sea-bosons is given by( q 7^ 0 ), 

4+q/2Ck-q/2 = (■^) ^ [Ak(q)ak(-q) + ak(q)^k(-q)] 

+Ti{k, q) al+q/2_qi/2(qi)ak-qi/2(qi - q) 

qi 

- T2{k, q) Y «Lq/2+qi/2(qi)«k+qi/2(qi “ q) ( 7 ) 

qi 

here, 

Ti (k, q) = i/l - hk+q/2 a/ 1 - hk-q/2 (8) 

72 (k, q) = A/Uk+q/2fik-q/2 ( 9 ) 

Ak(q) = Y^hk+q/2(l - hk-q/2) 

here, the sea-boson commutes with the total number of fermions, 

K(q), At] = 0 

and also the operator ak( 0 ) = 0 . Further, 

nk = + II«Lq/ 2 (q)«k-q/ 2 (q) - 5]aUq/2(^)“k+q/2(q) 

' ' q q 

and 

n^(k) = _ - _ 

exp{P{ey: - p)) + l 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 
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Also fik = {n-k)- The expectation value is with respect to the full interacting ground state. This quantity 
depends on the interactions that are present in the system and must be evaluated self-consistently. In fact, 
there is a deeper reason for introducing this. The exact formula for the number conserving product of two 
fermi helds and the sea-bosons may be expected to involve the number operator itself under the square root 
sign. This is made exceedingly likely by analogy with the bose case, where the square root of the number 
operator in the zero momentum state appears. In Appendix B the manner in which this exact correspondence 
may be deduced is hinted at. At this stage, it is pertinent to merely write down a formula for the sea-boson 
annhilation operator in the RPA-limit. The sea-boson is defined analogous to the condensate displacement 
boson, except the fermi case is more complicated due to the presence of the fermi surface. The sea-boson 
may be defined as follows (the rest of the details including a ’’proof’ of this fact and how it fits into the 
fermi-bilinear-sea-boson correspondence is relegated to Appendix B), 


ak(q) = 


(N) ’ 


(14) 


here 0(k, q) is a c-number phase that serves to randomly cancel out troublesome terms : this is also related 
to the ’’random phase” of the Random Phase Approximation of Bohm and Pines. Thus the above formula 
for the sea-boson is in the ’’Random Phase” approximation. 

This correspondence recovers the salient features of the finite and zero temperature aspects of the free 
theory provided we make the following assumption, the sea-bosons do participate in the thermodynamic 
averaging but come with an infinite negative chemical potential. This means that as far as the free theory 
is concerned, the average value of the sea-boson occupation is zero in the non-interacting case. The kinetic 
energy operator in the sea-boson language is given by, 


K = ^— al{q)aki(i) + N eo (15) 

k,q 

here eo is the kinetic energy per particle. Therefore, 


(alc(q)ak(q)) = 


1 


= 0 


(16) 


expiPi^ - pb)) - 1 

here —pB = oo. However, when there are interactions in the system, the answer is likely to be different. In 
particular, it is likely to be a non-analytic function of the interaction in such a way that it vanishes as the 
coupling goes to zero(this is demonstrated explicitly in Appendix D). Roughly speaking we may write. 


(aj,(q)ak(q)) « {^Pexp{-l/v) 




(17) 


where v is the coulomb repulsion parameter and V is the volume. All these do come out naturally from the 
correspondece written down above as we shall soon see. We have thus written down a useful correspondence 
between fermi and bose operators that recovers the salient features of the free theory at zero and finite 
temperature and it is clear that this correspondence is all that is needed to write down model hamiltonians 
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with any sort of interactions, like coulombic, with phonons e.t.c and extract exact non-perturbative (more 
precisely, non-analytic in the couplings) solutions. These solutions possess features that are impossible to 
capture via diagrammatic means let alone mean-field theory. Thus a strong case is to be made for this 
method as a new paradigm for condensed matter physics. 


III. FIELD OPERATOR IN TERMS OF DENSITY AND ITS CANONICAL CONJUGATE 


In this section, we introduce the canonical conjugate of the fermi/bose density distribution. The reason 
for doing this is that we would like to express the field operator itself in terms of the density and its 
canonical conjugate and consequently in terms of current and densities. None of these ideas are really new. 
For example, Sunakawa et.al.0 use the velocity operator as a canonical conjugate of the density in their 
investigation of the properties of He-II. The velocity operator is somewhat related to the current operator 
but is not exactly equal to it. The reason is that the current operator behaves like the conjugate of the 
density as far as commutation rules with the latter is concerned, but does not commute with members of its 
own kind(it is difficult to say this in words but will soon become clear). Let us postulate the existence of 
the object n(x(T) as the canonical conjugate of the density, 

[n(xCT), p(ycr')] =i S{x- y)S^y (18) 


[n(xcr),n(yCT )] = 0 


(19) 


It is clear that redefinitions of this object by amounts that involve translations by (more or less arbitrary) 
functionals of the density are not going to spoil the nature of the commutation rules above. However, we 
shall take the point of view that H is defined to be that(almost unique) object that satisfies the relation 
below (making mathematically rigorous sense out of all this requires the use of functional analysis and will 
be attempted in Appendix E). 


p(xcr) 
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i5n(xcr) 


( 20 ) 


Observe that p{x(t) = '01(xct)^(x(t) (technical problems involving the multiplication of operator-valued 
distributions at the same point are alleviated by assuming that we have the whole system in a box, with 
periodic boundary conditions on the fields, making any infinities only as large as the volume of the box itself, 
please refer to Appendix E for more details). Observe that (valid for both bose as well as fermi systems). 


[p{xa),i;{x'a')] = -(5^(x - xV<,,„'V'(xcr) 


( 21 ) 


rewriting this as a differential equation, 


( 22 ) 
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This may be solved (exponentiation of commutation rules is the more technical term) as, 


^/>(xcr) = ea;p(—in(x(T))F([p]; xcr) (23) 

Observe now that p = ijj. Therefore, 

^'''([p];xo-)F([p];xcr) = p(xcr) (24) 

this may in turn be solved and the final Density Phase Variable Ansatz (DPVA for short) may be written 
as, 

^/>(x(t) = e-* ‘^(['’l'’^'")(p(xa))5 (25) 

It may be noted above that redefinitions of 11 consistent with it being the canonical conjugate to p may be 
absorbed by a suitable redefinition of the phase functional $. Therefore, Eq.( |^) is in fact quite general. 
The crucial point of this whole exercise is that the phase functional 4) determines the statistics of the field 
'0(xcr). It may be shown (the proof is rather tedious and since this issue is not central to the practical 
computations, we defer the proof for a future communication) that imposition of bose/fermi commutation 
rules on ij) involves imposing the following restriction on the form of 4>. 

^([{p(yiCTi) - <5(yi - x')(5^iy}];xcr) 

+4>([p];x'cr') - $([p];xcr) 

- ^([{p(yiO’i) - ^(yi - X a) = rmr (26) 

where m is an odd integer for fermions even for bosons. This recursion is to be satisfied for all (xti) ^ (x ct ). 
It will be shown later that the restriction is far more severe, brought about by the need to recover the free 
case properly. It may puzzle the reader that the above statement implies that a random choice of the phase 
functional that ensures that the recursion is satisfied does not sufhce. This is mysterious, but is clarified 
by a conjecture in Appendix E. This is done by relating the canonical conjugate to the current operator 
and rewriting the DPVA in terms of current and densities. Again this type of idea has been addressed in 
the paper by Goldin et. al.0. However, many in this field(pardon the pun!) continue to be under the 
mistaken impression that the formula for the annhilation operator (say the fermi operator) in terms of the 
corresponding currents and densities depends on whether the fields in question are free or whether there 
are interactions in the system. This is shown to be false in the bose case, by demonstrating that there is 
a unique $ namely <I> = 0 that reproduces the free theory properly. Interactions just change the form of 
the hamiltonian but do not affect the form of the field operator in terms of currents and densities. The 
same is true but not easily seen in the fermi case; indeed throught this article we find that the bose case 
is much simpler and we shall take refuge under this rigorously justifiable edifice when confronted by fermi 


systems. Further, the formulas for the field operators suggested by Goldin, Menikoff and Sharp in their 
famous papeiS are according to our results only partially correct, since they have not actually introduced 
the phase functional d) and computed it (this will again become clear soon). 

Let us now write down a formula for the current operator in terms of the canonical conjugate and density, 

J = (^)(V'^V^-(V#V') (27) 

using the DPVA Eq.( psj) , 

J(xcr) =p(V$)-p(Vn+[-i$,Vn]) (28) 


From this it possible to deduce a formula for the conjugate in terms of currents and densities, 

;xcr)-y dl [-i$, Vn](y(T) (29) 

The line integral is along an arbitrary path from a remote point where all quantities may be set equal to 
zero. The field operator may now be rewritten exclusively in terms of currents and densities like so, 

here is canonically conjugate to the total number of fermions/bosons([Xoo-, = i S^^') = 

Sk''k(T''ktT The need for this is clear. The gradient of 11 does not involve the object Xq„ , when in fact 
it should. To put it differently, the field operator when commuted with should produce itself, whereas if 
we omit the object Xqo- then we find that the field operator commutes with the total number, which should 
not happen. These nuances are not very important for the practical computations as we shall see. It will be 
shown later that for bosons $ = 0 is the only possible choice and for fermions $ has to be fixed by making 
contact with the free theory. Uniqueness is assumed for the fermi case by making analogy with the bose case 
for which uniqueness may be proved. 


n(xcr) = Xq^ J (-i/p(yo'))J(yo') + 


IV. MAKING CONTACT WITH THE FREE THEORY 

In this section, we write down the kinetic energy operator in terms of the sea-displacements and determine 
the undetermined phase functional $ in the fermi case. The reason why the phase functional $ = 0 in the 
bose case will also be addressed here. Let us take the bose case first. It is clear at the outset that the 
choice $ = 0 satisfies the recursion Eq.( for bosons when one assumes that m = 0, an even integer. 
That this is the only possible choice is not at all clear. In order to verify this, let us write down the kinetic 
energy operator in terms of the density and its conjugate and show that an expansion in terms of the density 
fluctuations recovers the correct form of the dynamical density correlation function of the free theory(just 
the bose case). 
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( 31 ) 


K = [ [p(Vn)^ + ] + c — number 

J 2m 4/9 ‘ 

It may now be verified that an expansion in terms of density fluctuations leads to a hamiltonian that describes 
free harmonic oscillators, that may be easily diagonalised, it may also be shown that this diagonalised form 
reproduces the correct dynamical density correlation functions. The expanded form of the operator in fourier 
space is reproduced below for convenience, 


if = + (32) 

q^O 

A different choice of $ does not reproduce the free theory correctly. This is attested to by a simple calculation 
made in ID. Let us assume a form. 


^([p];^;) = 27r 



dy 9(x 


y){p{y) - po) 


(33) 


where 6{x) is the Heaviside step function. The above form clearly satisfies the recursion but does not 
reproduce the free theory as may be easily verified by the reader. 

The fermi case is somewhat more difficult. The difficulty is due to the fact that we must have a choice 
of 4) yf 0 that satisifes the recursion at the same time reproducing the free case. We shall take the point of 
view that the simplest choice for <i> namely linear in p should suffice. In any event, for the scheme to have 
practical significance, it is important for $ to be a simple functional of the density. We fix the coefficient 
in this ansatz by making contact with the free theory. Let us focus on the case of spinless fermions. In 
what follows, we restrict ourselves to zero temperature and a weakly non-ideal system, in this case, we are 
allowed to replace the Uk in the definition of Ak(q) by its noninteracting value at zero temperature. More 
interesting situations arise when the quantity fik is evaluated self-consistently, but we shall relegate these 
issues to future publicationsil. From equation Eq.( it is clear that redefinitions of the phase functional 
by amounts that do not depend on the density for example, $([p];x) ^ $([/ 9 ];x) -|- /(x), do not affect the 
formula for the field operator in Eq.( |^). Therefore, let us try the following ansatz for 4). 


®([p];x) = (34) 

qv^o 

Let us now write down the kinetic energy operator for fermions using the results of the first section. The 
kinetic energy operator was written down as, 

K ^'^—al{q)ak{q) + N eo (35) 

k,q 

It has been demonstrated in Appendix C that if one uses the form of the density fluctuation operator obtained 
by dropping quadratic terms in the sea displacements (the existence of such quadratic terms are hinted in 
Appendix B), this reproduces the RPA dielectric function. Since we know from prior experience that RPA 
is exact in the ultra-high density limit, we can use these two pieces of information to deduce a formula for 
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C/q(x) in terms of the properties of the free theory. First let us write down the RPA-form of the density 
fluctuation operator, 


= Xl[^k(q)ak(-q) + Au(-q)0|{^(q)] 


(36) 


where, 


Ak(q) = iynk+^/2(i^-^dZ^^ 


and the corresponding conjugate variable may be written down (that is, If in fourier space), 

1 


^q=(- 


2iN e, 


-)X][^k(-q)wk(q)ak(q) - Ak(q)u;k(-q)at(-q)] 


(37) 


(38) 


the dispersion is given by i.Ck(q) = From this the fermi field operator may be written down as, 


^(x) = e-* ^Ax)gi c^2(x)^ 


(39) 


where. 


C/i(x) = ^e* ‘i-lq 

q#0 


(40) 


^2(x) = ^ C/q(x)pq (41) 

q#0 

Using these facts, let us compute the equal-time version of the propagator below in the bose language and 
in the usual fermi language and equate the two expressions. In the sea-displacement language it comes out 
as. 


,t)) 


Po e 




(42) 


where. 


5'k.q(x) = -e * 'i-^(^-^^)Ak(-q)wk(q) + i [/q(x)Ak(-q) = -/k,q(x) (43) 

In the original fermi language it is, 

(V^t(x,t)^(x ,t)) = l^e*q-(-'-)0(fc;- |q|) (44) 

q 

Set Uq(x) = e“* 'i’^[/o(q) and I7o(q) is real. In order to derive a formula for Uo{q), let us equate the 
logarithm of the two expressions. 


log{{ipHx,t)'ilj{x.' ,t))) = log{po) + 


EK 

k,q#0 


1 

2 Af eq 


)2(]iq)2 + ([^^(q))2](A^(_q))2(g.q.(x-x') 

m 


1 ) 
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= logipo) + log{l + 1 - |q|)) « logipo) + 1 - l)B{kf - |q|) (45) 

q#0 


This leads to the following formula for the coefficient, 




wi(q) = ( 


4iVe2 ^ 


E(^)’(A^(-<!))= 


(46) 

(47) 


W2(q) = (^) ^(Ak(-q))^ (48) 

k 

In fact, in principle we could go all the way back to the expression in Eq.( ^0[ ) and say that we now have a 
unique correspondence between the fermi field operator and the corresponding currents and densities. In the 
next section, we write down and diagonalise the hamiltonian of interacting systems. It is shown that when 
only the lowest order sea-displacement terms/condensate displacement terms are included, it amounts to 
doing RPA/Bogoliubov theory. This hamiltonian is diagonalised in the fermi and bose cases and the single¬ 
particle spectral functions are computed. The bose case comes out nicely since, it is just the bogoliubov 
theory, but in the fermi case, we have to take extra care in properly diagonalising the hamiltonian in order 
not to lose the particle hole mode, the collective mode being more obvious. 


V. SPECTRAL FUNCTION OF INTERACTING SYSTEMS 

Let us make the following observation for future reference, 

RPA/Bogoliubov ^ Leave out the Quadratic Part in Eq.( 126) and Eq.( |3|) (49) 

It is pertinent at this stage to remark on the physical meaning of the above relation. In case of bosons, 
it is simple to visualise. Bogoliubov’s theory is exact provided there are large number bosons in the zero 
momentum state so that we may legitimately replace the number operator by its c-number expectation 
value. Also it is important that the system be weakly interacting so that the fluctuations of the number 
operator in the zero momentum state are small compared with its macroscopic expectation value. In the 
fermi case an analogous statement would be that the momentum distribution be sufficiently different from 
zero or unity for all values of the momenta. Also the fluctuations of the momentum distribution must be 
small. Thus for the fermi system, our approach gives good answers even for strong interactions that drive 
the momentum distribution away from zero or unity for all momenta so long as the fluctuations around 
these nonideal averages are small. In any event, the philosophy is that we have an exactly solvable class 
of models that describe correlation effects in many different contexts and this alone merits attention and 
serious study. In the end experiments may have to be used to ’calibrate’ these models so that they become 
a true description of the low-energy real world. 
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Part A 


Let us focus on the Lose case first. The bogoliubov hamiltonian may be written down by following the 
prescription of Eq.( ^ ). 

Hbog = ^ ^ [\/^f^-q/2 (“Q) + C^q /2 (q) \/^] [\/^^q/2 (q) + C^lq /2 (“^l) \/^] 

k q#0 

(50) 


In the above equation Nq is an operator, therefore this is the non-local bogoliubov hamiltonian. But we shall 
assume that it is legitimate to replace it with its c-number expectation value. It would be interesting to see 
what corrections to bogoliubov theory come about by incorporating this square-root of the operator. These 
correction terms tell us that fluctuations of the number of particles in the condensate are important and 
lead to correlations beyond the bogoliubov theory. This is in addition to correlations coming from quadratic 
terms that the prescription Eq.( neglects. When these approximations are implemented, and a further 
approximation Nq « A is made, it becomes exactly the bogoliubov theory introduced by Bogoliubov and 
Bogliubov and ZubarevS. It may be diagonalised quite easily, 


Hbog — l^q/q/q 
q 

and. 


(51) 


Z UJa 


2 U!, 


+ P0Wq,l ,+ , . 

-^)^dLq/2(-q) 


(52) 


dq/ 2 (q) = ( ^q + ^q + ^°^q 
Z 0J„ 



( 


-^q + Cq + PO'Cq 

2 Wq 



The dispersion is given by, 


(53) 


t^q = ^/fq + 2pol’qeq (54) 

here, po is the density of bosons in the condensate(not the overall density). From this one may deduce the 
filling fraction and dynamical structure factor, 


FILLING FRACTION 


/O = Nq/N = 1 - (1/A) ^(dj^/2)q(q)'^(l/2)q(q)) (55) 

q 

in other words, 


/o = Nq/N = 1 - (l/27rV) 



j 2/~^q + *=q + Po^^q, 

dq q (--^) 

Z UJa 


(56) 
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here p is the total density of bosons including those that are not in the condensate. 

DYNAMICAL STRUCTURE FACTOR 


S>{qt) = (pq(t)p_q(0)) = fVo([d_(i/ 2 )q(-q)(t) + d(i/2)q(q)W]M(l/2)q(q)(0) + rfL(l/ 2 )q(-q) (0)]) (57) 

in other words. 


^^(q, t) = iVo( —) exp{-i Wqt) 


(58) 


This method is truly powerful when applied to compute single-particle properties. The single-particle green 
function is difficult to obtain using conventional digrammatic methods or otherwise (See Kadanoff and Baym 
Ref.( ® ). For this one must first write down the field operator in terms of the condensate-displacements. 


n(x) « (^-^) ^ ea;p(fq.x)[dq/ 2 (q) - d(_q/ 2 (-q)] 


and the expression for the field operator is, 


V’(x) 


p-i n(x) 


V~p 


(59) 


(60) 


The propagator (all propagators in this article are evaluated at zero temperature, this means we may set 
the chemical potential equal to zero in the bose case) may now be computed and shown to be equal to the 
free propagator at ultra-high density. The interacting case is more interesting. The time-evolved version is. 


V’(x,t) 


^-i n(x,t) 


Vp 


and 


n(x,t) = (^-^)^ea;p(fq.x)(Aq-l-Rq)[/qe 


(61) 


(62) 


(V’^(0,0)V'(x,t)) = p ( e* n(o.o)g-* n(x,i) 


(63) 


In order to ensure that the free case is properly recovered, we use this somewhat illegal trick, but a trick 
that should be very palatable to most physicists, namely multiply and divide by the free propagator and in 
the division use the free propagator predicted by the bosonized theory and in the numerator use the free 
propagator obtained from elementary considerations. 


(V'’^(0,0)V'(x,t)) =exp [ fr 


where, 


AM 


. _ + Cq + PO'^q 'll 

^ 2 ^ 


(64) 


( 65 ) 
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Similarly, 


,-a;q + eq + po?^q,i 

^■' = <-2TZ- 


_ n(x,t) i n(o,o) 


(V’(x,t)'0f(o^O)) = p(e 


( 66 ) 


= exp [ ](r/'(x,t)r/>i'(0,0))/r 


/q(x, i) = (e-* - l)(^q + Bq)2 - (e"* - 1) 


here, 


From Kadanoff and BaymEI the spectral function may be deduced as follows, 

THE SPECTRAL FUNCTION 


(67) 


( 68 ) 


A(p,a;) = / dx 


^ + 00 -1 -1 


and, 




k.Xg-Z Ekt 


(69) 


(70) 


The above answer is the exact answer for the spectral function provided bogoliubov’s theory is adequate. 
Now we move on to the fermi case which is far more interesting and important. 

Part B 


In order to compute the full propagator for these systems, it is desirable to first ascertain, under what 
conditions these formulas are going to be valid. The answer is given by the assertion in Eq.( ^). Thus these 
answers for the single-particle properties are valid in the same limit in which RPA/Bogoliubov’s theory is 
exact. The assertion in the bose case in Eq.( has been verified. In order to verify the analogous assertion 
in the fermi case, we have to diagonalise the full hamiltonian given below (The fact that the RPA dielectric 
function comes out naturally from the prescription in Eq.( will be demsonstrated in Appendix C). In the 
fermi case, we have to diagonalise the full hamiltonian given below, 


H = ^Wk(q)aj^(q)ak(q) + X![^k(q)ak(-q) + Ak(-q)aj^(q)][Ak'(-q)ak' (q) + A,,'(q)aj^, (-q) 

q#o k.k' 


k,q 


(71) 


here Wk(q) = ^~Ak(—q). The zero temperature case is somewhat special, here we may assume that the 
sea-boson annhilates the non-interacting fermi sea, this means that we have to introduce a factor of Ak(—q) 
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in the dispersion that makes the kinetic energy operator positive definite. In order to diagonalise this we 
proceed as follows. Assume that the diagonalised form looks like this, 

i? = ^J}i(q)6^(q)&,(q) (72) 

i,q 

here 6i(q) and t>J(q) are ”dressed-sea-displacement” operators. The objects i take on values from an index 
set. The size of this set is the big issue here. Is it finite or does it have the same size as the number of 
points in k-space, or is it equal to the number of points on the fermi surface ? We shall find that answers 
to these questions are hard, and may be addressed only after coming to an agreement as to what sort of 
physics we hope to capture. Indeed, in many cases in physics one is forced to bend the rules or reinterpret 
mathematical formulae in order to capture what one is looking for. We find that we have to resort to such 
methods here as well. In particular, we find the following general feature, (i)i(q) are the roots of the RPA- 
dielectric function. Now the RPA-dielectric function is a complex quantity, as it is usually introduced in the 
textbooks. Therefore finding roots cannot mean finding the zeros of both the real and imaginary parts at 
the same time for this gives no root, both the real and imaginary part cannot be zero simultaneously. This 
leaves us with the following options, reinterpret the zeros of the RPA-dielectric function to be the maxima of 
the dynamical structure factor, in which case one gets both the particle-hole mode as well as the collective 
mode. The better option is to delay taking the thermodynamic limit until after all the summation over 
momenta have been performed. Then assume that the density is high enough and at the very end go to the 
thermodynamic limit, this ensures that both the particle-hole mode and the collective mode are properly 
recovered. These are admitedly difficult issues to grapple with, and the authors have attempted a different 
approach to deal with them. However, the traditional viewpoint on this matter is presented in the paper by 
Castro-Neto and Fradkinlli. The diagonalisation proceeds as follows, 

4(q)]ak(q) - ak(-q)]aj,(-q) (73) 

k k 

the corresponding inverted formula reads as, 


ak(q) =^[ak(q),6l(q)]&*(q)-^[ak(q),6i(-q)]6j(-q) (74) 

i i 

The quantities [6i(q), ak(—q)] and [ak(q), ^i(q)] are c-numbers and real. The i here could span a continuum 
(particle-hole mode) or be finite( actually there is just one, collective mode ). The diagonalisation continues 
unabated, 


[^*(q),ajc(q)] = ( 


Ak(-q) 


)ff*(q) = [ak(q),&I(q)] 


(75) 


[6*(q),ak(-q)] 


Ak(q) 


)£'*(q) 


(76) 
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( 77 ) 


[ak(q),fci(-q)] = ( , ^kq )g»(~q) 

w^(-q) + 


= [y Mk-q/2)-nHk + q/2) . 

The eigenvalues (I’i(q) are zeros of the real part of the RPA dielectric function. The RPA dielectric function 
is written down below, 


/ 1 , Wq Y^ni.(k + q/2)-ni.(k-q/2) 

= 1 + y y-- 

k m 


(79) 


As it stands, the above sum is ill-defined. In particular, if one takes the thermodynamic limit at the outset, 
and treats the above sum as the principal part, then one gets the real part of the RPA dielectric function. 
On the other hand, if one defers the taking of the thermodynamic limit until the very end, and instead takes 
the high density limit first, then one obtains the particle-hole mode as the argument below will attest. Let 
us rewrite the sum in the RPA dielectric function as. 


eflPA(q,0) 



A^(q)-A2(-q) 


1 _l_ y A^(q) 

V y (I>-f Wk(-q) 


V 


E 

k^ki 


A|(-q) 

0 - a;k(q) 


_ A^.(-q) 

V w - wki (q) 

Let us now assume that the volume V is fixed and we now go to the high density limit (fc_F 
equivalently when , |q| << k/), then we find, due to the fact below, 

Ak(—q) = 0 ; unless |k| « kf and k.q > 0 


(80) 
cx), or 


(81) 


The total number of terms in the above two sums is a small fraction of the total volume and as kp keeps 
increasing, the fraction gets smaller and smaller until it becomes small compared to unity and may be 
neglected, this means. 


. A^^(-q) 

^V^Uiiq)- Wki (q) 


= 0 


(82) 


from this we may deduce that particle-hole mode as. 


c^*(q)=^k.(q) + (^)A^.(-q) (83) 

as is clear from the above derivation, two points must be borne in mind, one is, we have to defer the taking 
of the thermodynamic limit until the very end, the other is to exploit the property of the object Ak(q), 
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namely, if |q| << kf, and Ak(— q) = 1 (Ak(— q) = 0,1 always) then |k| kj. Alternatively, we can solve 
for (I’i(q) as shown below. 




^k.(q) + (^) 


AL(-q) 


i + ^E 


A^(q) 


V ojk, (ql+i^kl-q) 


_V 

-al V ^ 


k^ki 


A£(-q) 

(q)-k;k(q) 


(84) 


We shall find these formulas useful later on when we try to write down the propagator. The collective mode 
in ID and 3D may be written down as shown below. 


( ■. /Ilk {kf + q/2)^-{kf-q/2)^exp{-X{q)) 
u!c-iD(q) = (— )\ - 


m 


1 - exp{-X{q)) 


(85) 


Kq) 


m Vn 


( 86 ) 


we may also write. 


kfq^2 ,,2,0, 


uJc-ioiq) = i^) + eg + 2 eq{^)coth{ 

fit lit Zi 


In 3D it is more familiarl(only for coulomb repulsion). 


i^c-3D(q) — tt>p[l + 


3 {.qvff 

10 


(87) 


( 88 ) 


For more general forms of interaction in 3D the answer may be obtained by computing the roots of the 
equation below, 


1 - 


nQ{vcfl^)lm 


1 .3 


{1 +-w[-r((7U_F) - Eq]} — 0 


L0‘- ' 5 

In 2D, the answer is not available in the books and may be deduced after some algebra as, 

-fc.f|qh/i I 27/// 


t^c-2D(q) = 


/^mi + HIIOl) 


ATz/m I /27r/m\2 

After all this, it is relatively simple to deduce the full propagator. For reference the free propagator is 


(89) 


(90) 


(V’^(x,t)V^(x',k)) = 

(?(>(x',k)V'^(x,t)) = poe”^wo/kA-')/k,<,(x') 

X ^k«7)0-O 

/k.q(x) = e*'i-^(^-^^)Ak(-q)wk(q) +i C/-q(x)Ak(-q) 


(91) 


(92) 

(93) 
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5k.q(x) 


(94) 


= _g-»q.x^ 


^2 N 


-)Ak(-q)wk(q) +i C/q(x)Ak(-q) = -/k,q(x) 


and, 


Zn = e 




^(Ak(-q))^(‘^k(q))^gi J2 


(95) 


The time evolved field operator in the interacting case is, 

'ip\x,t) = exp{ 17^^q(x)6j(q)e* “Aq)t)ea;p(_ ^ [/k;q(x)&,(q)e"* (96) 

k,q^0,i k,q^0,i 


where. 


t^k.q = /k,q(x)[ak(q),fel(q)] +/k,-q(x)[ak(-q),6,(q)] 


(97) 


7^o = exp{- Y /k,qW/k,q(x)[6»(q),o[(q)][ak(q),6l(q)]) 

k,q,i 

xexpi-^ Y /k,qW/k,-q(x)[ak(-q),6*(q)][ak(q),6l(q)]) 


X exp{— — Y /k,q(x)/k,-q(x)[ak(-q),5,(q)][ak(q),6l(q)]) 


k,q ,2 


The two full fermi propagators may be written down as, 
(V^t(x,t)^(x',t')) = 


(^(x ,t )'tp^x,t)) = |7^o^Zoppoe^‘‘•‘^ 


2,^ 


(98) 


(99) 

( 100 ) 


Again, it is desirable to use the trick we used in the bose case, namely multiply and divide by the free 
propagator and in the division use the form predicted by the bosonized theory and in the multiplication, use 
the form predicted by elementary considerations. This procedure also ensures that inspite of the fact we have 
not verified that the fermi fields written down in terms of the bose fields anticommute, the anticommutation 
rules are forced on the propagators by the free propagators which we know anticommute in the right fashion. 
This leads to the following forms for the propagators, 


,t)) = |7^oP|-Hol'‘e 


E 


)e* 


i(q)(t -t) 




(V'^(x,t)V’(x',t'))/r 


( 101 ) 


{tp^x = |7^o|■ 


l^ol^e^x 


: )C/k.q(x)e* 




Ek 


, /k.q(x)/k,q(x )e' 
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(V'(x ,t )i;^X,t))free 

( 102 ) 



In the above formula, the index i runs over both the collective mode as well as the particle-hole mode(z = c, k^) 
The momentum distribution may be evaluated in a different way by computing the expectation value of the 
number operator in Eq.( |l^ . This leads to the following answer. It includes contribution from both particle- 
hole mode and the collective mode. In Appendix D, we show how to derive the same momentum distribution 
using the equation of motion approach(actually just the collective part, for purposes of illustration). The 
full momentum distribution including the particle-hole mode is given below. 

(4ck) = 6{kf - |k|)Ei(k) + (1 - 6{kf - |k|))E2(k) (103) 


Fi{k) = 1- 


1 —nF(k-bqy 2 


l,q 


(ci,(-q) + ^^ + eq) 


9r (-q) 




^^2(k) = ^ 


nF(k-q) 2 


2 ,q 


(wi(-q) 


k.q 


- Co 


9^ (-q) 


(104) 


(105) 


In the above sum over i, one must include both the collective mode and the particle-hole mode(i = c,k^). 
A more general result is possible for systems that are significantly more nonideal. This comes about when 
one does not use the zero-temperature non-interacting values in the fermi-bilinear sea-boson correspondence. 
The form of the momentum distribution suggested by this is given in Appendix D. It is now very easy to 
write down a criterion for the breakdown of fermi-liquid behaviour. It is given by equating the step at the 
fermi surface to zero(the quasi-particle residue). 


Zf=F,{kf)-F2ikf)=0 (106) 

In the end, it is pertinent to address the claim made in the abstract namely that we are able to capture 
short-wavelength behaviour. The real issue here is that we have two length scales, one is the inverse of 
the fermi momentum the other is the lattice spacing. When one speaks of short wavelengths, one means 
wavelengths comparable to the lattice spacing. In the ultra-high density limit, where all the answers we 
have been deriving are valid, the inverse of the fermi momentum is much too small (compared to the lattice 
spacing) for the wavelength of any external field to be comparable to it. In other words, even if you have 
an external field that varies so rapidly in space that it changes sign every lattice spacing, the effective field 
induced by such an external field is still described by the RPA. To put it yet another way, the RPA is 
exact in the ultra-high density limit. Some have argued that this limit is uninteresting since in this limit, 
the coulomb interaction is completely screened out and therefore in this regime we just have a fermi liquid. 
We find that this argument is not entirely true. In fact, we have shownEl that when the inverse of the 
fermi momentum is small compared to the lattice spacing, it is still possible to increase the value of the 
dimensionless coupling strength sufficiently so that fermi liquid behaviour is destroyed. Actually, our results 
turn upside down almost all conventional wisdom, first we find that fermi liquid behaviour persists in ID 
for sufficiently weak coupling strengths (even when we assume the interactions are hard core delta-function 
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interactions), in contrast to the Lieb-Mattis solution of the Tomonaga-Luttinger model. Second, we find 
that fermi-liquid behaviour breaks down in more than one dimension for sufficiently strong values of the 
coupling strength in contrast to the answers obtained by Castro-Neto and FradkinEl. In fact we find that 
fermi liquid behaviour persists in all three dimensions for sufficiently small values of the coupling strength 
and is destroyed in all three dimensions for sufficiently large values of the coupling strength. It may be 
argued by the reader that our results are not foolproof either, for one, we have neglected several terms in 
the hamiltonian and those terms are small only in the limit when RPA is exact. The other points are the 
technical shortcomings, such as the fact that we have not proved the fermi case as rigorously as the bose case, 
like the fermi commutation rules are not explictly verified e.t.c.. Notwithstanding all these shortcomings, a 
strong case is to be made for the revision of entrenched dogma about fermi and luttinger liquids. 


VI. CONCLUSIONS 

Let us summarise the results obtained so far. We have succeeded in reducing to quadratures the propaga¬ 
tors of both bose and fermi systems. We have also computed the momentum distribution of interacting fermi 
systems and written down a formula for the quasi-particle residue in terms of the electron-electron repulsion. 
From this we obtain a criterion for the breakdown of fermi-liquid behaviour. The results we obtain contra¬ 
dict widely held views about ID systems, in particular the Lieb-Mattis solutioni of the Tomonaga-Luttinger 
model suggests that the momentum distribution of a ID system with delta-function interactions exhibits no 
discontinuity at the fermi momentum. This is in contrast with the exact formula above that does in fact 
exhibit such a discontinuity for sufficiently weak values of the coupling strength and is destroyed only for 
larger values of the coupling strength. We attribute this discrepency to flaws in the linearised dispersion 
model (i.e. Tomonaga-Luttinger model) and the ensuing algebraic manipulations, not all of it transparent. 
In any event, the authors feel that the Tomonaga-Luttinger model is not a caricature of the homogeneous 
interacting fermi system in the ultra-high density limit. Luttinger Liquid theory is based on the unproven 
assumption that the low-energy behaviour of the homogeneous interacting fermi system in one-dimension 
is correctly described by the exactly solvable Tomonaga-Luttinger model. Our results show that this is 
false. The important qualitative features of the homogeneous interacting fermi system namely the presence 
or absence of a fermi surface cannot be surmised by examining the properties of the Tomonaga-Luttinger 
model, especially when the interactions are weak. 

APPENDIX A 

In this appendix we prove some assersions made earlier. First the definition of the condensate displacement 
annhilation operator. 

c^q/2(q) = i^=)blbq (107) 
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In order to define the quantity O = in a manner acceptable to most physicists, we proceed as follows. 

O is defined to be that operator that commutes with the number operator Nq 


[0,No] = 0 


(108) 


in the basis in which Nq is diagonal and possesses non-zero eigenvalues (not an unreasonable assumption 
considering the fact that even in the most strongly interacting systems Nq is macroscopic), call them {IVg}, 
then the matrix elements of O in the same basis are going to be 1/ ^/Nq. Having thus provided all the 
matrix elements, the definition of O is complete. We have to now show that (iq/ 2 (q) satisfies canonical bose 
commutation rules. The simplest way of doing this is to use the polar decomposition of bo 

bo = exp{-i Xo)^/l^ (109) 

here Xq is the hermitian operator canonically conjugate to Nq = ^o^o, that is, [Xq, Nq] = i. This decompo¬ 
sition correctly reproduces the bose commutation rules of bo and 6 q. For example, 

[6o, &o] = bobl - b^bo = exp{-i Xq) Nq exp{i Xq) - Nq = 1 (110) 

This means that (iq/ 2 (q) = z^bq, where Zq = exp{i Xq). Since, [zoj^q] = 0 and [^Oj^q] = 0, and [zoXo] = 0, 
it follows that dq/ 2 (q) and 6q both satisfy the same commutation rules since, Zg now behaves effectively as 
a c-number (as regards commutation rules with 5q, fc^^and zq. It is worthwhile pointing out this fact, 

Mq/2(q),^o] 0 


rather, 


[dq/2(q),iV] = 0 (111) 

though not obviously so. In order to prove this, 

K/2(q),^^] = K/2(q),^^o] + K/2(q), ^q'Vl 

qVo 

= [exp{i Xg), Xg]6q -f eXp{i Xg) ^ [6q, 6q,&q'] 

qVo 

= {exp{i Xg)Xg - NoeXp{i Xg))6q -I- expii Xg)6q 


= [i Xq, No]exp{i Xo)bq + exp{i Xo)bq = —exp{i Xg)6q -I- exp{i Xg)6q = 0 


( 112 ) 


Next, one would like to prove the equation Eq.( For this we simply plug in the definition Eq.( IW) into 
Eq.( ® and verify that is reduces to an identity. The details are as follows, first define, 
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— -^0^k,0^q,0 + [<^k+q/2,o(\/^)c^k( —q) + <^k-q/2,0C^k (q) ( V^^)] 

+ ^(l/2)(k+q/2)(^ + q/2)d(l/2)(k-q/2)(k “ q/2) (113) 

The proof involves these cases, 

(i) k = 0 and q = 0 
In this case, 

Lq^o = Nq = 6q6o (114) 


(ii) k + q/2 = 0 but k — q/2 ^ 0 


l^k= —q/2,q ('s/l^o)<^—q/2( q) — 


(115) 


(iii) k — q/2 = 0 but k + q/2 ^ 0 


l^k=q/2,q *^q/2 (*l) ( 's/^0 ) 5q6o 


(116) 


(iii) k — q/2 ^ 0 and k + q/2 ^ 0 


ik.q = 'l(l/2)(k+q/2)(l^ + q/2)d(l/2)(k-q/2) (k - q/2) = bl^^^^eXp{-i Xo)eXp{i Xo)5k-q/2 = ^k+q/ 2 ^k-q /2 

(117) 


Therefore in all cases, 

-bk.q = ^k+q/2^k-q/2 (US) 

and thus Eq.( |^) follows. Finally, we would like to clarify the finite temperature case. In particular, what 
is the chemical potential of the condensate displacement bosons ? Is it zero or is it the same as that of 
the parent bosons ? The answer may be found by computing the thermodynamic expectation value of the 
number of bosons in the condensate A^O) 

{No) = N- Xmi/2)q(q)'^(l/2)q(q)) (US) 

q#o 

We also know the answer from elementary considerations, it is. 
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1 


( 120 ) 


{No)=N-Y, 




exp{P{eq - n)) - I 


where p is the chemical potential of the parent bosons. Then it follows that, 


1 


(rf(i/ 2 )q(q)'^(i/ 2 )q(q)) - exp(/3(eq - m)) ” 1 


( 121 ) 


In other words, the chemical potential of the condensate displacement bosons is the same as that of the 
parent bosons. 


^parent f^cond/displ 


( 122 ) 


APPENDIX B 

In this Appendix, we try to to make plausible the correspondence between the number-conserving product 
of two fermi fields and the sea-bosons. Let us rewrite the bose case (Eq.( I)) more suggestively, 

^k+q/2^k—q/2 — O(^) '^q.O + [•\/^k+q/22lk(~Ct) + ^k(q)'\/^k—q/2] 

+ X! "^k+q/2-qi/2(qi)^k-qi/2(-q + qi) 


In the bose case, 


qi 


^ ^k-q/2+qi/2(qi)^k+qi/2( q + qi ) 


qi 


(123) 


^k(q) — ^k-q/2,0'^q/2(q) 


(124) 


and. 


0(k) = N 4,0 


(125) 


Observe that the suggestively extravagant notation in Eq.( 123) is meant to imply that a very similar relation 
holds in the fermi case which we reproduce below, 


4+q/2‘^k-q/2 — 0(k) (5q,0 + [V’^k-|-q/2^k( —q) + ^k(q)V^k-q/ 2 ] 


+ X! ^k+q/2-qi/2(qi)^k-qi/2(-q + qi) 

qi 

~ ^k-q/2+qi/2(qi)^k+qi/2(~q + qi) (126) 

qi 

Here Ak(q) depends on two momentum labels unlike in the bose case. This has to do with the fact the 
now 0(k) no longer has the simple structure we saw in the bose case. We must now invert this relation 
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and obtain a formula for the operator Ak(q). It is not at all clear that this object will behave like an exact 
boson annhilation operator. The alternative is to write down an ansatz for an exact boson in analogy with 
the bose case and determine the unknown in the formula by imposing canonical bose commutation rules. 

Qk(q) = -4-q/2^(^^>q)<^k+q/2 ( 127 ) 

y^k-q/2 

The unknown operator M (k, q) has to be related to some number-conserving fermi bilinear by demanding 
that the operator ak(q) obey canonical bose commutation rules, 

[ak(q),ak'(q')] = 0 (128) 

[ak(q),aj^/(q')] = (5k^k''^q,q' (129) 

It is at present beyond the authors to arrive at a formula for M(k, q). Notwithstanding this, it is still useful 
to capture some sort of an approximate correspondence like the one introduced in Sec. II. The relations 
written down there have the following positive features 

(i) They recover the RPA dielectric function at zero and finite temperatures. 

(ii) They capture the correct four-point and six-point functions at zero and finite temperatures. 

(iii) The formula for the sea-boson in Eq.( 0) when plugged into the correspondence for the number operator 
in Eq.( ^ gives an identity. 

The only negative aspect of this correspondence is that 

(I) The mutual commutation rules between the off-diagonal fermi bilinears is recovered correctly only upto 
terms linear in the sea-bosons. That is, somehow the operators on the right side of these commutations rules 
should not be too different from their approximations obtained by dropping terms higher than the linear 
order. This is no doubt a strong assumption. This is in fact equivalent to RPA(perhaps even better than 
RPA). 

The definition of the seaboson is incomplete without the prescription of the phase 0(k, q). In order to 
derive an expression for this, we again make heavy use of the bose case which we have proved rigorously in 
Appendix A. There we found that plugging in the expression for the condensate-displacement boson into the 
correspondence resulted in an identity in case of q 0 (the q = 0 case being special). This identity comes 
about in a very specific fashion. In the general form of the correspondence outlined in Eq.( |123[) , we find 
that the sum on the right that comes with a negative sign is identically zero( for q ^ 0 ) and the sum on 
the right that comes with a positive sign is equal to the left-hand side, except in ’’rare” cases when either 
k + q/2 = 0 or k — q/2 = 0. We shall adopt the same approach in the fermi casse and try to fix the phase 
0(k, q) such that the identity is satisfied in the manner just described. Let us now write down the potential 
identity, 
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4«, ‘"•-■‘W-,,. 


+Awf-q)c^ C-- g(k,q)( »^(J^-q/2) a„ i_ 

+ ^k( qjC^+q/^e ( ) Ck_q/2^_^ 


+7’i(k,q)4+q/2( ^ 

qi#q,o 


v-^ n 


.^(k + q/2 qi) ^; e(k-qi/2,-q+qi)„-i 6l(k+q/2-qi/2,qi)\ 

(N) > 


Ck-q/2 


/.X 1 If /''T-^(k + q/2). 1 ,n^(k — q/2), 1 

-T2(k,qK-q/2 \.n )"( /^n 


' V"k-q/2 V’^k+q/2 {N) 


(TV) 


X ^ nk_q/2+q,e*"(‘^+‘i^/2’-‘i+'i^)e-* 

qi/q,o 


0(k-q/2+qi/2,qi) 


( 130 ) 


Here, since we are not involved in proving the rigorous correspondence, but just the salient features, we are 
entitled to some leeway. In particular, we shall turn a blind eye to the fact that there exist these objects 
Ti(k, q) and T 2 (k, q), in fact set them both equal to unity, just for the moment. The exact correspondence 
in terms of the Hk(q) seems to suggest exactly this. Then we find that, if we choose our d(k, q) to be such 
that 


0{k - qi/2, -q + qi) = 0{k + q/2 - qi/2, qi) 


( 131 ) 


and. 


E 


Tik-q/ 2 +qie* »(k+qi/2.-q+qi)g-* 0 (k-q/ 2 +qi/ 2 .qi) ^ q 


qi#0,q 


( 132 ) 


then all is well. Terms that were linear in the sea-bosons are vanishingly small in the thermodynamic limit, 
and are important only when both the sums on the right side are identically zero for some reason, that is, it 
is ’’rarely” important just like in the bose case. It is not really important to write down an explicit formula 
for the phase function 0(k, q), it is merely sufficient to show that it does what is required of it, namely, it 
provides the ’’random phase” that cancels terms that enable the whole machinery to run smoothly. Lastly, 
we have not yet verified that this sea-boson obeys canonical commutation rules. This is again a tricky 
problem, it is likely to be resolved by the exact approach which is beyond the scope of this article. It is 
merely sufficient to point out that this is likely to come about due to the strong likelyhood that the phase 
d(k, q) is actually a functional of the number operator. 

The correspondence that we have just defended is nothing but a more elegant version of the correspondence 
introduced by the pioneers like Castro-Neto and FradkinEi. Any criticism that may be leveled against our 
approach may equally well be leveled against theirs. The only difference between our approach and theirs 
is that the single-particle properties which they are so fervently seeking are far more elegantly recovered by 


26 














our approach since we do not linearise the bare fermion dispersion or use the clumsy Luther const ructionlliJ. 
Indeed, we have even shown that the answers for the ID case are different from the Tomonaga-Luttinger 
model that linearise the bare fermion dispersion. 

The other issue worth addressing at this stage is the validity of the prescription in Eq.( It can be 
seen from the exact correspondence in Eq. ( |l26| ) that as q ^ 0 terms that correspond to corrections to the 
RPA-form of the full hamiltonian vanish at least as fast as |q|/fc/. The RPA-terms themselves don’t vanish 
and tend toward (limq_o Ak(—q) ^ 0 as in the bose case) 

^ Vnir-^k(-q)(133) 

k k 

In order for the prescription in Eq. ( |^) to be accurate, it is important for the interaction Uq to possess these 
properties, first it must vanish for large enough q(or small inter-particle separation), 

Lim|q|^co kf Vq^O (134) 

where cq is small compared to unity and positive. This ensures that the only possible contributions come 
from small q where corrections to the RPA-form themselves are small. In addition, if we also make sure that 
the interaction vanishes fast enough for large inter-particle seperations so that, 

LzTO|q|^0 Vq ^ |q|^ (135) 

where D = 0,1, 2, ..(larger the better), then our formalism is in fact EXACT as kf oo(or sufficiently 
large). It may be argued that this state of affairs is most likely uninteresting since it may not be realisable 
in practice, when it is, it merely leads to a fermi liquid. This is a valid point. But it is worth pointing out 
that non-fermi liquid behaviour can still emerge in such systems when the interaction strength (with the 
same functional form) becomes strong enough. These considerations also tell us that for an interaction of 
the delta-function type in ID, provided the strength is weak enough, we have a fermi-liquid in contrast to 
the Lieb-Mattis solution of the Tomonaga-Luttinger model. 

In any event, the philosophy is, that having introduced sea-bosons, we more or less forget about the fact 
that it was fermions that motivated their introduction in the first place and instead try to write down a whole 
new set of models in terms of the sea-bosons calibrate them appropriately so that they capture the salient 
features of the real world. It is not a tautology to remark that we have in our hands a whole class of exactly 
solvable models of correlated fermions that is easier to use than mean-field theory itself but capture effects 
significantly beyond diagrammatic perturbation theory, like the nonanalytic dependence of the momentum 
distribution on the coupling strength(written down in Appendix D). 

APPENDIX C 

In this appendix we demonstrate that the RPA dielectric function is recovered exactly by selectively 
retaining parts of the coulomb interaction that lead to RPA. We know that the kinetic energy in the bose 
language is given by, 
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(136) 


For this let us choose, 


where, 




kin = ^(^)a[(q)ak(q) 

k,q 




q#0 


2V 


PqP- 


(137) 


P<i = E[^k(q)ak(-q) + Ak(-q)a[(q)] (138) 

k 

From this it may be shown that the RPA dielectric function is recovered as the following demonstration 
shows. Assume that a weak time-varying external perturbation is applied as shown below, 

= ^(C/e,i(q,t) + u:,,(-q,t))pq (139) 

q#0 


here, 


C/e,,i(f,t) = [/oe*‘i-"-“‘ 

Let us now write down the equations of motion for the various bose fields, 

*^(ak(q)) = Wk(q)(aL(q)) + (y )Ak(-q) (-q)(aL'(q)) + V(q)(ak'(-q))] 

k' 


(140) 


+ (C^exi(q, i) +t^e*rt(-q.0)Ak(-q) 


(141) 


-*^(«k(-q)) = Wk(-q)(ak(-q)) + (y )Ak(q) ^[A^/(-q)(a^, (q)) -1-Aj^/(q)(a^t (-q))] 


k' 

+ (Uext(qjt) -j- U^^f( — q, t))Ak(q) (142) 

Now, let us decompose the expectation values as follows, 

(ak(q)) = Uext(q,t)Ck(q) + U*^t(-q,t)Di,(q) (143) 

(aLt(_q)) = u:,^(-q,t)C^(-q) + C/,,,(q, t)74* (-q) (144) 

The coefficients C'k(q) and D^(—q) satisfy, 

w C'k(q) = a;k(q)C'k(q) + (^)Ak(-q) Y^[K' (-q)C'k' (q) + V (q)A'k' (“q)] + Ak(-q) (145) 

k' 


-w L>k(-q) = Wk(-q)L»k(-q) + (y)Ak(q)Xl[V(q)-Dk'(-q) + Ak'(-q)C'k'(q)] -f Ak(q) (146) 
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Now, the effective potential may be written as, 


t) = t) + {y){p^q}' Uextifl, t) (147) 

here, 

(P-q) = C^ext(q,t)(P-q)' + (-Q, t) (p_q)" (148) 


Using the fact that, 


(P-q)' = Xl^k(-q)C'k(q) +^Ak(q)£)^(-q) 

k k 


(149) 


Solving these equations and using the fact that the dielectric function is just the ratio of the external divided 
by the effective potential we get, 


_ Uext{^.,t) _ 1 , ;((q nF(k + q/2) - nF(k - q/2) 
C^e//(q,t) 


(150) 


Which is nothing but the RPA dielectric function of Bohm and Pines. 


APPENDIX D 


In this appendix we use the equation of motion approach to solve for the momentum distribution and 
compare it with the solution obtained via exact diagonalisation as described in the main text. The equations 
of motion for the bose propagators read as. 


(*^ - Wk(q)) 


-i(TaJ,(q)aJ^,(q )) 


(Tl) 


= 5, 


k.k'^ 


q,q 


m 


+ (-^)Ak(-q)^[Ak-(-q)--+ Ak"(q)-^riT^^-^ 

k" 


(151) 


d . ,,-i(ra[‘(-q)a^(q')) 


(*^ + Wk(-q)) 


(Tl) 


= -(y)Ak(q)Xl[V'(q) 


-i(Taj^(,(-q)aj^,(q')) -t(Ta^„ (q)aj^, (q')) 


(Tl) 


+ ^k" (“q)‘ 


(Tl) 


(152) 


The boundary conditions on these propagators may be written down as, 

-t(Ta|{^‘(-q)aL(q')) -i(Ta|{,^‘”*^^(-q)aL (q')) 


(Tl) 


(Tl) 


(153) 


-i(TaJ,(q)aL(q')) -i(Ta^i^ "^\q)al(q)} 


w 


(Tl) 


(154) 
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(155) 


Sit) = i—^) ^M^nt) 




-i (3' 

The boundary conditions imply that we may write, 
-t(ra^(q)aL (q')) 


(Tl) 


= X! ea:p(w„t) 


-r(Tai((q)at,(q')) 

(Tl) 


(156) 


(157) 


-i(Taj^*(-q)aj^, (q')) 
(Tl) 


exp{uJnt) 

n 


-i{Ta^^{-q)al, (q')) 
(Tl) 


(158) 


and, uJn = (2 TT n)l(3. Thus, 


[iUn - Wk(q)) 


-t(Ta"(q)a|^, (q )) _ (5i, i,o5q 


(Tl) 


—i [3 


+ (y)Ak(-q)2][Ak"(-q)--+Ak"(q)- 


(Tl) 


(Tl) 


(159) 


(iw„ + Wk(-q)) 


-i(Ta|,”(-q)aj^,(q')) 

(Tl) 


= -(y)Ak(q)y3[Ak"(q) 


-i(Taj^”(-q)aj^, (q')) 
(Tl) 


+ Ai, 


'(-q)-^ 


(q) 


4'(q')), 


(Tl) 


(160) 


Define, 


X] Ak(-q) 


-r(TaO(q)a;^,(q')) 

(Tl) 


Gi(q, k',q';n) 


(161) 


X] Ak(q) 


-i(Taj^”(-q)aj^,(q')) 

(Tl) 


G 2 (q, k',q';n) 


(162) 


Multiplying the above equations with Ak(—q) and summing over k one arrives at simple formulas for Gi 
and G 2 . 


Gi(q,k ,q ;n) = Ak'(-q) 


q.q 


-i (3iiujn - Wk'(q)) 


+ fn (q) [Gi (q, k', q'; n) + G 2 (q, k', q'; n)] 


(163) 


and. 
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G 2 (q,k',q';n) = /*(-q)[Gi(q, k', q'; n) + G 2 (q, k', q'; n)] 


(164) 


G 2 (q,k ,q ;n) = Gi (q, k , q ; n) 


Gi(q,k ,q ;n) + G 2 (q,k ,q ;n) = Gi(q,k ,q ;n)/(l - /*(-q)) 


. 1 . (l-/«(-q))Ak'(-q)<5q,q' 

Gi(q,k ,q ;n) = (—^)- 


(1 - /^(-q) - /n(q))(*w„ - Wk'(q)) 


(165) 


G 2 (q,k',q';n) = (^-;) 


/«(-q)Ak'(-q)^q,q' 


-* /3 (1 - /^(-q) - /n(q))(*w„ - Wk' (q)) 


(166) 


Gi(q, k , q ; n) + G 2 (q, k , q ; n) = (—^) 


V(-q)'^q,q' 


-*/? (l-/^(-q)-/«(q))(*w„-a;k'(q)) 


-*(TaC(q)aJ^,(q')) 

(Tl) 


^k.k'^q.q' 


-z /3(zw„ - Wk(q)) 


1 ^^Uq^ Ak(-q) 


+ (^)(l7)7- 


Ak' (-q)<5, 


q,q 


'-i(3 V (zw„ - a;k(q)) (1 -/^(-q) -/n(q))(iw„ - Wk'(q)) 
The zero temperature correlation function of significance here is, 


(167) 


-i(al (q')ak(q)) 


(168) 


This may be obtained from the above formulas as 

- *(4' (q )ak(q)) = -(y )Ak(-q)Ai^/ (-q)'^q,q' 


duj 


C 2tt i {ioj - a;k(q))(*w - (q))(l - /*(-q) - /„(q)) 

(169) 


where G is the positively oriented contour that encloses the upper half-plane( upper half-plane, because we 
need (aj^, (q )ak(q)) and not (ak(q)a|^, (q )) ) Thus the problem now reduces to computing all the zeros of 
(1 ~ /n(~q) ~ /n(q)) fkat have positive imaginary parts. It may be shown quite easily that. 


eflPA(q, iw„) = 1 - /*(-q) - /„(q) 


In ID, the dielectric function is evaluated as follows, 

1 m ikf + q/2)^ + {^? 

1 - /:(-.) - /.(.) = 1+= 0 


This leads to the root. 


m 


\ 


ikf + q/2Y - {kf - g/2)^exp(-(^)(^)) 

l-exp(-(^)(;k)) 


(170) 


(171) 


(172) 
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Therefore the final result may be written as, 


1 


h(q )ak(q)) = (-) 






here, 


V (ujR(q) + UJk(q))(uJit(q) + Wfc' (q))(^)(^)2(f)2ujR(q)(cosh(A(q)) - 1) 


A(<,) = (^)(^) 

m Vq 


(173) 


(174) 


/A (kf+ q/‘^)^ - {kf - q/2)^exp{-X{q)) 

( m V- i-exp{-X{q)) - 


(175) 


In other words, 


(4cfe) = npik) + {2Trkf) j ^ 


Afe-gi/2(-9l) 


2LOR(qi){ixiR(qi) + uJk-q,/2{qi)?{^){cosh{X{qi)) - 1 ) 


/ + 00 
-oo 


dqi _ Afc+,^^/ 2 (-gi) _ 

27r 2ujR{qi){ujR{qi) + ^^k+q^/2^dl)Y{^){cosh{X{ql)) - 1) 


(176) 


Note that the above formula possesses a non-analytic dependence in the coupling strength {cosh{-^){^) — l) 
, an unmistakeable signature of a non-digrammatic result. Next, we would like to provide formulas for the 
momentum distribution when we use the correct interacting expectation values in the fermi-bilinear sea- 
boson correspondence. The results obtained from these formulas are likely to be very different from the 
weakly nonideal case, which in any case is not very interesting. The answers are given below, 

n^(k) S' 2 (k) 


Uk = 


^i(k) ^ ^i(k) 


here, 


5i(k) = l + ^ 


nk- 


q.* 


(w*(-q) 


-eq)^ 


'(-q) + E 


1 ^k+q 


(d'i(-q) + ^ + eq) 


:5*^(-q) 


(177) 


(178) 


^2(k) = ^ 


rik- 


(‘^i(-q) + ^-eq)^ 


:ff*^(-q) 


(179) 


also the form of the ”RPA” dielectric function and its zeros (I’i(q) are now different. The ”RPA” dielectric 
function is given by. 


^RPA\ 




(180) 


^^(q) = E 


^k-q/2 ~ ^k+q/2,_ 1 

Eq7i~i 


k (‘^A(q)-^)^ 

The commutators are given as before, except for three changes. In the new approach. 


(181) 
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(182) 


Ak(q) = i^nk+q/a)! - nk-q/ 2 ) 


Next, the zeros are slightly different. The collective mode has to be computed self-consistently, whereas the 
particle-hole mode may be written down as described earlier, 


w*(q) = Wk,(q) (-j^)--- 


AL(-q) 




_ Ik 

V 




V a;ki(q)+‘^k(-q) V Z^k^ki uJi^. (q)-u;k(q) 


(183) 


The last change is in the form of {7o(q), here we have to make sure we use the finite temperature noninteracting 
values. The other issue that is also of interest is whether the momentum distribution evaluated using the 
fermi-bilinear/sea-boson correspondence is the same as that suggested by the full propagator. We have found 
that the answer to this is difficult and probably in the negative. This does not mean that the whole program 
is wrong. Some comfort and confidence in these manipulations may still be retained by demonstrating that 
the expression for the number operator is consistent with the RPA-form of the fermi creation operator. 
Again here, we have to be content with a weak form of this requirement. We take the point of view that it is 
sufficient to show that the commutator between the total momentum of the electrons and the field operator 
comes out the same in both the original fermi language and in the sea-boson language. The total momentum 
of the electrons has the form, 


P=.^kcJ,Ck (184) 

k 

In the sea-boson language, it takes the form, 

P = X! ^ a[(q)ak(q) (185) 

k,q 

Therefore in the original fermi language we have. 


[P,i/'^(x)] = f V 


(186) 


In the sea-boson language we have, 

[P,V’^(x)] = ^q[oJ,(q)ak(q),V’^(x)] 

k,q 

= Xl^®k(^)(“5k,q(x))V'^(x) + ^q(-6f£ q(x))V’^(x)ak(q) = i V V'^(x) 

k.q k,q 

as it should be. All this points to the fact that the answers for the momentum distribution and propagators 
should not be taken too literally, rather one must be content with the qualitative predictions that are most 
likely accurate, that also seem to contradict conventional wisdom. 

APPENDIX E 
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In the late 70’s and early 80’s, attempts were made to write down field theories that describe scalar mesons 
in terms of observables like currents and densities rather than the creation and annhilation operators. The 
motivation for doing this stems from the fact a theory cast directly in terms of observables was more 
physically intuitive than the more traditional approach based on raising and lowering operators on the Fock 
space. This attempt however, raised a number of technical questions, among them was how to make sense 
of the various identities connecting say the kinetic energy density to the currents and particle densities and 
so on. Elaborate mathematical machinery was erected by the authors who started this line of research^ 
to address these issues. However, it seems gaps still remain especially with regard to the crucial question 
of how one goes about writing down a formula for the annhilation operator (fermi or bose) alone in terms 
of bilinears like currents and densities. The bilinears in question namely currents and densities satisfy a 
closed algebra known as the current algebr ai. This algebra is insensitive to the nature of the statistics of 
the underlying fields. On the other hand, if one desires information about single-particle properties, it is 
necessary to relate the annhilation operator (whose commutation rules determine the statistics) to bilinears 
like currents and densities. That such a correspondence is possible was demonstrated by Goldin, Menikoff 
and Sharpli. However, they have not explicitly written down such a formula nor have they clarified some 
important issues such as whether this formula changes when one consider interacting fields rather than free 
fields. The general beliefi is that these formulas are different for interacting fields. It is shown here that this 
is in fact not the case, interactions in the system merely cause a change in the hamiltonian but do not affect 
how the annhilation operator is related to local currents and densities. The attempts made here are partly 
based on the work of Goldin et.alEl Ligouri and Mintchev on generalised statisticsE^ and the series by Reed 
and SimouE. As has been demonstrated earlier, for the bose case we had to choose $ = 0. We argued that 
this choice was unique. In the fermi case the choice was different but was also unique due to the necessity 
of recovering the free theory. In this section, we write down a mathematically rigorous statement of this 
uniqueness criterion. This exercise also settles the issue regarding the delicate question of multiplying two 
operator-valued distributions at the same point and other related issues, like the meaning of the square-root 
of the density distribution. For this we prove this lemma. 

Lemma Let T be smooth a function from a bounded subset of the real line on to the set of reals. Also let / 
and g be smooth functions from some bounded subset of Vf' to reals. Let us further assume that the range 
of these functions are such that it is always possible to find compositions such as T o f and they will also be 
smooth functions with sufficiently big domains. They possess fourier transforms since they are well-behaved. 

If, 


•^( ./W ) = ff(x) 


(187) 


and. 


/(x) = ^/qe* 


(188) 
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(189) 


then the following also holds, 


k 


[^( ^/qr_q(k) )] 4.0 = 5k (190) 

q 

where T'q(k) = ea;p(q.Vk)- Here the operator Tq(k) acts on the k in the Kronecker delta on the extreme 
right, and every time it translates the k by an amount q. 

Proof 

Proof is by brute force expansion. We know, 


Hy) = E 


jF(")(0) 


n—0 


nl 


(191) 


therefore, 


Ti /(x) ) = .F(O) + E 

n—1 


jP(»)(0) 

n\ 


E /qi/q2--/qnea:p(i(Eqi)-x) 

{q.} *=i 


= Ee*'^"5k (192) 

k 

This means (take inverse fourier transform), 

^(0)4,. + f; ^ 

"=i {qil 

= ffk (193) 

This may also be cleverly rewritten as 

^(0)4,0 + E -E ./qx/q.-/q„r_q,(k)T_q,(k)...T_q„(k)4,0 

"=1 { q .} 


= 4 (194) 

and therefore, 

4 = E /q^-q(k) )] 4,0 (195) 

q 

and the Proof is now complete. 

Now we would like to capture the notion of the fermi density operator. Physicists define it to be p{x) = 
'ilj*{x)'ijj{x). Multiplication of two fermi fields at the same point is a delicate issue and we would like to make 
more sense out of it. For this we have to set our single particle Hilbert Space: 
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n = Lp{n^)<^w 

Here, Lp{TZ^) is the space of all periodic functions with period L in each space direction. W is the spin space 
spanned by two vectors. An orthonormal basis for W 




A typical element of Ti is given by /(x) A basis for H is given by 


^ = hi ;^ea;p(iqn.x)(^^s;n = (m, 712 , 713 ) G Z^,s e {T,i}; 

We move on to the definition of the fermi-density distribution. The Hilbert Space " is the space of all 
n-particle wavefunctions with no symmetry restrictions. From this we may construct orthogonal subspaces 

= p+n®'^ 


O’ 


Tensors from TL 


1. The space 


= P_'H® 

O" 


^ are orthogonal to tensors from P 


The only exceptions are when 77 = 0 or 77 = 

/O’ 


is the space of bosonic-wavefunctions and the space P_ is the space of fermionic 
wavefunctions. The definition of the fermi density distribution proceeds as follows. Let v be written as 


The Fermi density distibution is an operator on the Fock Space, given a vector / O ^ ^ ^ single 

particle Hilbert Space, and a tensor ip in the n-particle subspace of of T{P), there exists a corresponding 
operator p{f O) v) that acts as follows: 


[p(/(^^^)‘/?]n(xi(Tl,X2Cr2, ....,X„Cr„) = 0 

ii p G P^ ” and 

n 

[P(/0^^)‘P]«(X1CT1 , X 2 (J 2 ,-, Xn(Jyj) — ^ ^ (^1^1 7 ^2^2 7 -: ^n*^n) 

2 = 1 

when p G P_ . The physical meaning of this abstract operator will become clear soon. Let us now define 
the current density in an analogous fashion. To Physicists, it is, 

J(x) = (^)[V’^(VV’) - (196) 

To Mathematicians it is an operator similar to the densityS. Given a typical element f associated with 
the underlying single-particle Hilbert space, there is an operator denoted by Jsif^v), ( s = 1,2,3 ) that 
acts on a typical tensor from the n-particle subspace of the full Fock space as follows. 
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[•^s(/(^t’)‘P]n(xi(Tl,X2(T2, ...,X„Cr„) = 0 


if tp G 7f®",and and, 


?;)V3]„(xiCTi,X2Cr2, ...,x„cr„) 


(197) 


= -i^{/(xfc)a(f7fc)Vg + -[Vj/(xfc)]a(a/e)}<^n(xiai,X2f72, ...,Xnan) (198) 

k^l 

if tp G 7i_ . For the bosonic current it is the other way around. Having done all this, we would now like to 

write the DPVA more rigorously. Now for some notation. As before, let g = exp{i km.x)0^r (the square 
root of the volume is not needed as we want all operators in momentum space to be dimensionless). Then 
as before. 


■0(k„r) = c{g) 


(199) 


p(k„r) = p{g) 


SpCk^r) = p{k^r) - 


( 200 ) 

( 201 ) 


js{kmr) =Js{g) 


( 202 ) 


Sjsikmr) = js(k„r) 


Having done this, we would like to write down another formula for the canonical conjugate. 

Vn(x(T) = (—1 /p(x(t))J(x(t) + V<I)([p]; xcr) — [—i $, VH] 


(203) 


(204) 


Then we have(bear in mind here that we have distinguished between the c-number Ad? and the operator 
p(0r) whose expectation value is At)?). 


(t Hm) — (^^ 0 ) 1 , 1 


At)? 1 + ^ ^P(knf)7k„(qm) 


] ^qm.o + F([p];q„r) (205) 


where. 


exp{i qm.x)F([p]; q^r) = V$ - [-i $, VH] (206) 

As regards the object Aor that is conjugate to the total number is concerned, we must retain it as it is, 
since, it will ensure that the total number when commuting with the field operator is the field operator itself 
rather than the incorrect answer zero. For q^ 7 ^ 0 
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1 




1 


Ek„'5Kk«»’)Tk„(qm) 


^qr 

Pri 


i-'^j(Pn*")7p„ (qm) ] l^qm,0 


i qm.F([p];q^r) 


(207) 


In order to define Xor in terms of fermi fields, we have to make use of the fact that this object does not 
commute with the total number of fermions. This means it cannot be expressed exclusively in terms of 
number-conserving fermi bilinears like currents and densities. This will mean that we merely invert the 
formula in Eq. ( and solve for Xor as, 


Xor = - + * Y ^^[(\/^ ■i/i(qnr)r_q^(km))(iV° + y^^p(qnc)Tq„(k^)) 2exp{-i Y 

km q„ q„ q„ 


X (/.([p];q„r)rq„(k„))] 4^.0 (208) 

Define an operator which is defined to be the formal expansion that the formula itself suggests, 

'0(k„r) = exp{-i n)Xq^r) exp{i 

qm qm qm 

(209) 

We would now like to write to write down a statement that would require a proof. This conjecture when 
proven will vindicate the DPVA. 

Conjecture 

There exists a unique functional $([p];xr) and a unique odd (for fermions, even for bosons) integer m such 
that the following recursion holds, 

^([{p(yic^i) - ^(yi - xV<ti,a}1;xct) 


+^([p];x'cr') - $([p];xcr) 


- ®([{p(yicri) - <5(yi - x)(5^i,„}];x cr ) = mn (210) 

and has the following additional effects. The domain of definition of '0(k„r) ( in which the series expansion 
converges ), is the same as that of ^(k„r) and it acts the same way too. In other words, 

tpiknr) = ipiknr) ( 211 ) 

We know how the ingredients of V'(k„r) namely the current j(k„r) and the density 5p(q„r) act on typical 
elements of the Fock space, and we know how '0(k„r) acts on the Fock space, we just have to show that 
the complicated ^(k„r) acts the same way as the simple ip(knr). Moreover, this is true for a unique phase 
functional d). Lastly, we would like to defend the above “fourier gymnastics” by pointing out that the real 
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space formulation is not well-defined due to the fact that the line-integral that appears in the formulas is 
difficult to define, any attempt is equivalent to the above approach. The other reason for attempting a 
rigorous formulation is the fact well-known to mathematicians that it is not possible to have a self-adjoint 
canonical conjugate of a positive definite self-adjoint operator. Since p is positive definite, the natural 
question that arises is whether 11 self-adjoint ? We take the naive physicist’s approach to this issue, namely 
we allow for sign changes in p and argue that these merely amount to translating the phase functional <I> by 
constant amounts, thus not altering the general framework. Within this framework. If is indeed self-adjoint 
and all is well. It is also worth remarking that the overall conjugate 11 has two contributions, one from the 
position independent part Xq^, and the other is from terms involving currents and densities. The latter 
contribution is manifestly self-adjoint. The possible lack of self-adjointness of the overall conjugate stems 
from the canonical conjugate of the total number which cannot be expressed in terms of fermi bilinears. 
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